1994 


N95- 19012 




— i 

V»*j 


f) 


J 

r 


6 


NASA/ASEE SUMMER FACULTY FELLOWSHIP PROGRAM 


MARSHALL SPACE FLIGHT CENTER 
THE UNIVERSITY OF ALABAMA 


NUMERICAL FORMULATION OF COMPOSITION SEGREGATION AT CURVED 
SOLID-LIQUID INTERFACE DURING STEADY STATE SOLIDIFICATION 

PROCESS 


Prepared By: Jai-Ching Wang, Ph.D. 

Academic Rank: Associate Professor 


Institution and 

Department: Alabama A&M University 

Department of Physics 


NASA/MSFC 

Office: 

Division: 

Branch: 


Space Science Laboratory 

Microgravity Science & Application Division/ES 
71 

Electronic & Photonic Materials Branch/ES 75 


MSFC Colleague(s): Sandor L. Lehoczky, Ph D. 

Frank Szofran, Ph.D 
Dale Watring 


XLV 




INTRODUCTION 


The lateral solute segregation that results from a curved solid-liquid 
interface shape during steady state unidirectional solidification of a binary 
alloy system has been studied both analytically and numerically by 
Coriell, Boisvert, Rehm, Sekerka (1). The system under their study is a two 
dimensional rectangular system. However, most real growth systems are 
cylindrical systems. Thus, in a previous study(2) we have followed Coriell 
etc. formalism and obtained analytical results for lateral solute segregation 
for an azimuthal symmetric cylindrical binary melt system during steady 
state solidification process. The solid-liquid interface shape is expressed as 
a series combination of Bessel functions. In this study a computer 
progi 'll-', has been developed to simulate this lateral solm segregation. 


FORMALISM 


In this section we present the basic equation and boundary condition 
used in this calculation. The diffusion equation for an azimuthal system is 


n 3 2 c'(r‘,z') [ i 3c'(r',z') ,3 2 c'(f,z') v 3c'(f.z;) n 
3r‘ 2 <" 3r' ~d~ z u 


— ( 1 ) 


Where D is diffusivity of solute in the liquid, V is the velocity of 
solidification, c',r', z',w' are dimensional solute concentration, radial and 
axial coordinates, and interface thickness. The boundary conditions are 


VC' I (k-l)=P(jS^2 - 

)>■, Sr- > (2) c(r,z) = c 0 . 

3C'(r',z') 


3r'i 


-=0 at r =R 


(3) 

(4) 


Where k is distribution coefficient. The variables can be 
nondimensionalized by letting c=c'/co, r=rVR, z=z'/R, w=w'/R and (3=VR/D 

where R is radius of the ampule. The diffusion equation and the boundary 
conditions become y 
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P c lfk Q- aC(r ’ z) dC(r,z)3w(r) 
^ ' dz\ 2"* 


dr\ dr 


c(r j 2 =oo)=i 


•(7) 


3C(r,z) n 

— - — 0 at r= 1 
dr 


6) 


( 8 ) 


In the limits of very small |3, by using variable separation method, 
equation (5) can be solved to give 

c(r > 2)=AJ 0 (ar)ef 1+[1+(2a/b>2l ’ ,Sz + A' —..(9) 

From boundary condition c(r,z=°o)=i^ We have A ,_.j 

dC(r,z)_ 


From boundary condition 

^HL-aJ l( ar)=0 
d r 


dr 


0 at r= 1 


, We have 


( 10 ) 


let u n be the zeros of Jl(°0, where n=l to infinite. 

2u 2 1/2 

[l+(~) ] 

By using p(n) to denote P . The general solution is 


°° R 1/2 

C(r,z)=l+l^€-P(z-wo) + I A n Jo(u n r)e7 tl+p{n)] z 
k n=1 


( 11 ) 


There is another boundary condition which needs to be satisfied, i.e. 

BCrfk l) =dC(*%z) aC(r,z)aw(r) 

^ 7 dz\ 


9ri dr ' 


( 12 ) 


The solid-liquid interface shape is assumed to be parabolic 5r 2 and 
is expressed as a series combination of Bessel function, i.e. 
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oo 


w(r)=w 0 + 2 5Cn)J 0 (u n r) 
n=1 


(13) 


dw(r) 

Assume both w(r) and dr to be small. We obtain the solute 
concentration at the interface to be 

pi(r 1 v 6(n)Jo(u n r) 




C sl (r,2)=kC| 

= 1-fVl V ^( n )Jp( u n r ) 

w n=i l+2k/[p(n)-l] (15) 

This shows that the solute concentration in the solid at the 
interface is a function of w(r), p and k. 

NUMERICAL CALCULATIONS AND RESULTS 
The results obtained from these calculations show that 

(1) The solute concentration at the interface at r= 0.7 which were 

calculated by using the general expression eq (ll).with n from 1 to i is 
denoted by j . The same value calculated by using the approximate 
expression eq(15).with n from 1 to i is denoted by The results in 

figure 1 show that if we have include more than 40 terms in the 
calculations, the two expressions give almost exact results. 

(2) .The solute concentration in the solid at the interface obtained 
analytically show that the compositional segregation in the solid is 
proportional to the deviation of the interface from planarity. The 
proportional factor being the product of P and (1-k). The solute 
concentration at the interface C si (x) calculated by using the general 
expression and by using the approximate expression for a small p limit 
agree very well.The results for p =0.173, k=4 and 5 = 0.05, 0.1, 0.2, 
0.375 and 0.4 are calculated and the result for 5 =0.375 is shown in 
figure 2.The results of C1 j were calculated by using the general 

expression eq (11) with n from 1 to 89, while the results of U2j were 
calculated by using the approximate 

XLV-3 



fr-MX 


1X91 " S U0 ^ * Wl 99s uopdB D X am*!* 



I. J0J ui |. 0jqBj uj jBqj sb j: rpu0d©p >f pue $ jbiiuits ©ABq s:qns©j 

9S9qx 'l=* PUB 0=X }B PHOS ©qj UI ( X )IS D ©OBJJ0JUI ©qj }B UOT5BJ5UOOUOO 

0^108 ©qj 0 ai 6 0 M ‘9 PUB ‘<3 ‘>f jo S 0 njBA snojJBA joj ‘j. 0 jqBj u| g 

‘snipuj sji jo uoipunj b sb ©dBqs 

9UIBS ©qj 9ABq n® S9pmnduir 10BJJ9JUI 3U0J9JJIP JOJ (x)! » 3 0OBJI0JUI 
9qj IB UOIJBJJU9DUOO ©jnios oqj pjqj Avoqs p ©jnSij ui sjins©j ©q X -p 

•p9J9pISU009J 9q 

o; p90u Z,£99£ 31 = ^ JOj sqnsoj ©qj j * $ ijbuis joj pqBA si qoiqM(£i)b© 
uoiss 0 jdx 0 9jBUiixoiddB ©q ? Suisn Xq P 9JB l nO l BO 9JB SJ I nS 9 J 9 S 9 q 1 
L£ 99£'3I P UB 6£99£3’I V£99£3I‘0= d -ioj ©objjsjui ©qj jo ©gp© 

©q5 l=x pus ©objj©jui ©qj jo j©ju©o ©q; ‘o=x qjoq jb uoijbia©p ©dBqs 
©objj©jui ©qj jo ©pniijduiB ©qj o } fBuopjodoid XjjBsuq si / x) ’is^ ©objj©jui 

oqj IB U0PBJ5U90U0D ©jrqos ©qj iBqj Moqs £ ©jngjj ui s^ns©j © qi *£ 

‘68 05 I uiojj u qjiM (gj)b© uoiss©idx© 








If 

P 

K ri 




P 













4 - 3.14 d(0) ri(1) 


If 

Jt 

p 

r* 



If 

p 

5 

n 

p 


0.1 

0.01 

0.05 0.9972 1.003 


Wk, 

4 * 3.14 

»(0) 

»( I ) 



4 * 3.14 


Si(I) 


0.1 

0.01 

0.1 0.9944 1.006 


2 

0.01 

0.05 1.0029 

0.997 


2 

0.1 

0.05 

l . Oi . 

0.979 


0.1 

0.01 

0.2 0.9887 1.011 


2 

0.01 

0.1 

1.0058 

■'>94 


2 

0.1 

0.1 

1.0326 

0.958 


0.1 

0.01 

0.3 0.9831 1.017 


2 

0.01 

0.2 

1.0117 

0.988 


2 

0.1 

0.2 

1.0652 

0.916 


0.1 

0.01 

0.4 0.9775 1.022 


2 

0.01 

0.3 

1.0175 

0.982 


2 

0.1 

0.3 

*8 

0.873 


0.1 

0.01 

0.5 0.9719 1.028 


2 

0.01 

0.4 

1.0233 

0.976 


2 

0.1 

0.4 

l . ljcr 45 

0.831 


0.1 

0.01 

0.6 0.9662 1.034 


2 

0.01 

0.5 

1.0292 

.97 



2 

0.1 

0.5 

1.163 

0.789 


0.1 

0.01 

0.7 0.9606 1.039 


2 

0.0 1 

0.6 

1.0350 

0.964 


2 

0.1 

0.6 

1.1957 

0.747 


0.1 

0.01 

0.8 0.95498 1.045 


2 

0.01 

0.7 

1.0408 

0.958 



2 

0.1 

0.7 

1.2283 

0.705 


0.1 

0.01 

0.9 0.9494 1.05 



2 

0.01 

0.8 

1.0467 

0.952 



2 

0.1 

0.8 

1.2609 

0.663 

r 

ft 1 

ftftl 

J ft ft Q437 1 ftS6 , 


2 

0.01 

0.9 

1 . 05 :: 

1946 



2 

0.1 

0.9 

1.2935 

0.62 







2 

0.01 

1.0 

1.0583 

0.94 




2 

0.1 

1.0 

1.3261 

0.578 



P 



|r 


p 


r 

c 



If 

P 

5 

c 

p 


Ir 

5 C d<0) C d(I) 


Wk 



o 

u si(0) 




K 

4 - 3.14 

«(0) 

1 ) 


4 - 3.14 



4*4.14 






0.1 
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0.01 
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0.1 0.9459 1.055 


10 
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Table 1. Solute segregation for solid-liquid interface shape w(r)= or 2 for 
various values of k, p, and S 
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